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ABSTRACT 

Galactic nuclei host central massive objects either in the form of supcrmassive black 
holes or nuclear stellar clusters. Recent investigations have shown that both compo- 
nents co-exist in at least a few galaxies. In this paper we explore the possibility of a 
connection between nuclear star clusters and black holes that establishes at the mo- 
ment of their formation. We here model the evolution of high redshift discs, hosted in 
dark matter halos with virial temperatures > 10^ K, whose gas has been polluted with 
metals just above the critical metallicity for fragmentation. A nuclear cluster forms as 
a result of a central starburst from gas inflowing from the unstable disc. The nuclear 
stellar cluster provides a suitable environment for the formation of a black hole seed, 
ensuing from runaway collisions among the most massive stars. Typical masses for the 
nuclear stellar clusters at the time of black hole formation {z ^ 10) are in the range 
10^ — 10^ Mq and have half mass radii < 0.5 pc. The black holes forming in these 
dense, high redshift clusters can have masses in the range ~ 300 — 2000 Mq. 

Key words: black hole physics - instabilities - stellar dynamics - galaxies: nuclei - 
galaxies: formation - galaxies: star clusters 
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1 INTRODUCTION 

Central massive objects are known to reside in the centers of 
a large fraction of galaxies either in the form of supcrmassive 
black holes (BHs) or nuclear stellar clusters (NCs). Both 
BHs and NCs have masses that correlate with those of their 
hosts (Ferrarese et al. 2006; Wenher et al. 2006), suggesting a 
symbiotic evolution between these nuclear components and 
the larger structure of their galaxy hosts. 

NCs, with masses, Mnc, ranging between 10^ — 10^ Mq 
(Walcher et al. 2006; Cote et al.(2006); Bokert et al. 2004; 
Balcells et al. 2007) appear to be more common in galaxies 
with stellar masses below ~ 10" Mq, while BHs (Mbh > 
10® Mq) preferentially inhabit galaxies more massive than 
10" Mq (Ferrarese et al. 2006; Wenher et al. 2006). De- 
spite this apparent dichotomy, a smooth transition exists 
between NCs and BHs. Some galaxies do host both a NC 
and an embedded BH (Seth et al. 2008; Graham et al. 2009) 
and observational estimates seems to indicate mass ratios 
Mbh /Mnc ~ 10"* - 10"^ (Seth et al. 2008). Interestingly, 
BHs coexisting with NCs have masses that covers the inter- 
mediate range (10'^ — 10® Mq) between stellar and super- 
massive BHs. 

Two main scenarios have been proposed to explain the 
origin of NCs. In the first, NCs are the end product of mul- 



tiple mergers of globular clusters following their sinking by 
dynamical friction into the central part of the host galaxy 
(Tremaine et al. 1975; Capuzzo-Dolcetta & Miocchi 2008; 
Lotz et al. 2001). In the second, NCs form in situ after 
inflow of gas from the outer part of the host (Milosavlje- 
vic 2001): this latter possibility seems more consistent with 
stellar population studies in NCs (Cote et al.(2006)). 

Whereas models of NCs in a cosmological context are 
still lacking, models of BH seed formation have been de- 
veloped in the ACDM cosmogony. How BH seeds form and 
how they evolve into the supermassive variety has been the 
subject of many studies (Begelman & Rees 1986). A nat- 
ural path to BH formation relies on the first generation of 
stars (Pop HI) born in metal free environments (Madau & 
Rees 2001; Volonteri et al. 2003). These stars are expected 
to be more massive than their current counterparts (Abel 
et al. 2000; Omukai & Palla 2003; Freese et al. 2008; locco 
et al. 2008). PopIII stars that form above a threshold mass 
of ~ 260 A'Iq are believed to directly implode leaving a seed 
BH of similar mass (Heger et al. 2003). Alternative routes to 
BH seed formation typically exploit gas-dynamical processes 
in metal-free galaxies (direct formation models) (Haehnelt 
& Rees 1993; Loeb & Rasio 1994; Eisenstein & Loeb 1995; 
Bromm & Loeb 2003; Koushiappas et al. 2004; Begelman, 
Volonteri & Rees 2006; Lodato & Natarajan 2006). Gravita- 
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tional instabilities can indeed lead to a vigorous gas inflow 
into the very central region of protogalactic discs, supplying 
the necessary matter for the formation of BH seeds. 

The above scenarios have a common fundamental in- 
gredient, i.e., the gas hosted in the progenitor halo has not 
been polluted by metals. Metal pollution marks the end of 
the formation of PopIII stars and of their seed black holes as 
fragmentation and formation of low mass stars start as soon 
as gas is polluted above a critical metallicity threshold, Zcth 
(Schneider et al. 2006; Bromm et al. 1999, 2002; Omukai et 
al. 2008; Clark et al. 2008, Santoro & ShuU 2006). PopIII 
stars with masses between 140 and 260 M0 are believed to 
be responsible for starting metal pollution, as they do not 
leave a relic black hole but explode releasing all metals in the 
surrounding medium. Star formation is also an impediment 
for strong gaseous inflows in the center of massive halos as 
it limits the available mass supply. This does not necessary 
means the end of the seed formation epoch as a new chan- 
nel for BH formation opens (Omukai et al. 2008; Devecchi 
& Volonteri 2009, hereafter DV09). 

DV09 showed that the first episode of low-mass, PopII 
star formation can provide the conditions for the forma- 
tion of NCs and BH seeds arising in the most compact of 
these clusters. Gravitationally unstable protogalactic discs 
whose gas is only mildly polluted by metals, can still pro- 
duce strong gas inflows into their central regions. Stars start 
to form in the central few pc where the gaseous mass has 
been accumulated and the densities are high enough for star 
formation to set in. Clusters formed in this way are crowded 
places. Star-star collisions in their core can proceed in a 
runaway fashion and a very massive star of ~ 1000 M© can 
build up before the flrst supernova explodes (Freitag et al. 
2006b,a; Portegies Zwart & McMillan 2002; Gurkan et al. 
2004). At low metallicities, the flnal fate of this very mas- 
sive star is to coUape into a BH with mass similar to its 
progenitor - as mass loss at low metallicity is negligible. 

One of the factors that regulate the competition, at a 
given redshift, between different formation scenarios is the 
metal enrichment history. Both observational and theoreti- 
cal works have been performed in order to investigate how 
metals are distributed as cosmic structures form and evolve 
(Scannapieco et al. 2003; Tornatore et al. 2006; Savaglio 
2006; Savaglio et al. 2005; Prochaska et al. 2003; Prochaska 
et al. 2007, Kulkarni et al. 2005; Li 2007). These studies 
have shown that metal enrichment proceeds in a very inho- 
mogeneous fashion. This opens the possibility that different 
modes of BH seed formation co-exist in different regions of 
the Universe at the same redshift, instead of being mutu- 
ally exclusive. As the mean metallicity increases with cosmic 
time, growing bubbles of polluted gas where the first NCs 
and their BH seeds form, co-exist with relatively pristine 
regions where PopHI star formation can still be feasible. 

This is the first of a series of papers in which we de- 
velop a model aimed at tracing self-consistently metal en- 
richment, following the hierarchical build-up of halos. We 
consider metal enrichment and radiative feedback from the 
first stars together with feedback of successive generations 
of stars. Our study is applied in the context of BHs and NCs 
formation in order to determine: (i) if one of the mechanisms 
outlined above is dominant in the overall seed formation 
scenario; (ii) if a clear transition in redshift exists between 
different formation channels or if they co-exist in the same 



cosmic epoch; (iii) the characteristic redshift at which each 
process reaches its maximum efficiency; (iv) when the first 
NCs are expected to form. 

In order to address the above questions, we develop in 
this paper the semi-analytical scheme, aimed at describing 
the evolution of a single halo in which a central massive ob- 
ject forms. We adopt a ACDM cosmology with fib = 0.041, 
Sim = 0.258, SIa = 0.742, h = 0.742 and = 0.963 as given 
by five-year WMAP data (Dunkley et al. 2009). 

The outline of the paper is as follows. In Section 2 we de- 
scibe our recipes for the formation and evolution of a gaseous 
disc, hosted in a high redshift dark matter halo. In Sec- 
tion 3 we set our model for the formation of NCs resulting 
from dynamically unstable discs and of BH seeds ensuing 
from runaway collapse of stars in NCs. In Section 4 we il- 
lustrate the results and compare the results of DV09 with 
those found with our new technique. Finally, in Section 5 we 
briefly summarise our model. 

In the subsequent papers we will apply our semi- 
analytical code to merger tree histories extracted from the 
code PINOCCHIO (Monaco et al. 2002a,b). This wiU allow 
us to trace halo evolution, keeping track also of their spatial 
positions, a condition necessary in order to understand the 
relevance of radiative and SNae feedback. 



2 DISC FORMATION 

We here consider a dark matter halo of mass Mh, virialising 
at redshift Zvir, assuming that it follows the density proflle 
of an isothermal sphere. Virial radius -Rvir, circular velocity 
Vii and virial temperature Tvir are inferred from the top hat 
collapse model (see for example Barkana & Loeb 2001). 

Gas falling into the potential well created by the dark 
matter is heated by adiabatic compression and shocks, at- 
taining a mean temperature equal to Tvir- We assign an ini- 
tial gas mass (fib/flm) Afh and assume that gas follows an 
isothermal density profile Pgas(r). 

It has been proposed that a critical metallicity ^crit 
exists above which star formation in the low mass mode 
starts. We here focus only on those halos whose gas has 
been enriched above Z^it- In the following we describe the 
details of our treatment for the evolution of halo gaseous 
component, star formation and BH seed formation. 

2.1 Gas cooling 

We assume the standard model for disc formation (White & 
Rees (1978), Cole et al. (2000), Mo, Mao & White 1998). To 
this purpose for a given Tvir and metallicity Z of the gas, 
the cooling timescale TcooI can be computed as 

_ 3 fimH ksTvir 
- 2 p,,,{r) A(r.i„Z)' 

where A(rvir, Z) is the cooling function. We adopt A(rvir, Z) 
tabulated in Sutherland & Dopita (1993). The cooling radius 
rcooi{t) at time t is defined by imposing Tcooi = t. Matter 
inside rcooi has time enough to cool down from its initial 
temperature. 

For the disc to form, gas has to both cool down and 
collapse. The timescale for disc formation will therefore be 
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the longest between the cooling and the free-fall timescales. 
In analogy with the computation of rcooi(t) we can define 
a free-fall radius rft (t) that sets the radius inside which the 
gas has sufficient time to free-fall into the central structure. 
During the initial phase of the collapse a mass of cold gas, 
A^coid, develops at a rate 



Mcold = 47rpgas(rmin(f))r,nin(i)'^n 



(2) 



where rmin = min[rcooi, ra]. 

The cold gas is assumed to condense into a rotationally 
supported disc whose angular momentum Jd is a fraction 
jd of the angular momentum Jh of the halo. If the specific 
angular momentum of baryons is conserved during collapse, 
id = Mcoid/Mh. 

The surface density profile of the disc follows a Mestel, 
isothermal profile: = T,q{Ro/R) where scale parame- 

ters Eo and Rq are calculated assuming mass and angular 
momentum conservation (see DV09]Q. 

In the following Sections we describe how scale parame- 
ters of the disc vary with cosmic time since (i) the disc mass 
keeps increasing due to the infall of fresh cold gas accret- 
ing from the halo (as described by equation [2]) and (ii) the 
gas content changes due to gravitational instabilities that 
triggers episodes of infiow and star formation (outlined in 
Section lO). 



The infiow rate that channels mass in the central region 
of the disc is then: 



2r]4 
G 

GV^ 



^ -1 



\iT^ g"^ qIYjqRq 



(5) 



Following DV09 we assume that the inner surface den- 
sity profile changes its scaling with radius as oc R~^^^ , ac- 
cording to simulations by Mineshige & Umemura (1999). 
The transition radius, i?tr, between the two profiles is given 
by: 



Rtv — 



47rEoi?o' 



(6) 



where Minf is the net mass infiowing from the outer disc 
{R > Rti, see DV09 and equation [9|. 

For a given set of Mh, Tgas and Qc, a maximum angular 
momentum exists below which a disc becomes unstable A < 
l/8Sli, /SlmQc x/Tvir/Tgaa where Tgas is temperature of the 
cold gaseous disc. Discs with higher rotation are more stable 
agains fragmentation as the amount of mass that needs to 
be removed from the Mestel disc in order to reach Q — Qc 
is lower. The inflow mass decreases for halos with large spin 
parameters. 



2.2 Gravitational instabilities 

The stability of the disc is described in terms of the Toomre 
parameter Q (Toomre 1964). As the surface density increases 
since more cold gaseous mass is added to the disc, Q de- 
creases and eventually drops below the critical value for disc 
stability Qc- We here assume Qc = 2 (see DV09, Lodato & 
Natarajan 2007, and Volonteri et al. 2008 for a discussion on 
Qc). At this point the disc develops bar-like structures that 
can lead to a redistribution of mass and angular momentum. 

Infiow of gas driven by gravitational instabilities can be 
described in terms of an effective viscosity i/. The infiow rate 
is given by (Lin & Pringle 1987, see also Lodato 2008): 



2ttuT,{R)\ 



dlnn{R) 



d\nR 



where Q is the angular velocity and 



V = -q 



TvGE{R) ■ 



(3) 



(4) 



We set in our reference model 77 = 0.3 (Gammie 2001) but 
we also compare our result with simulations run with 77 = 1 
(high efficiency inflow) and 77 = 0.1 (low efficiency infiow). 
The sound speed Cs corresponds to a gas with temperature 
of 8000 K0. 



^ Note that the evolution of the collapsing structure and the ini- 
tial assembly of the disc is described in a way similar to DV09. 
The main difference is that the assembly of the disc is here fol- 
lowed solving esplicitly for the time evolution of the cold gas, 
instead of assuming that a fixed fraction of the halo gas cools 
down instantaneously. 

^ In this paper wc consider only halos illuminated by a fiux in 



2.3 Star formation 

If the disc grows strongly unstable, star formation sets in, 
consuming part of the gas that would otherwise fiow into the 
central region. Further evolution of Minf is then regulated by 
the competition between star formation in the outer Mestel 
disc and tunneling of gas in its center. 

An unstable disc starts fragmenting into bound clumps 
once the gravitationally induced stress exceeds a critical 
value corresponding to Merit = 0.06 Cg /G (Lodato & Natara- 
jan 2007). We allow fragmentation in the Mestel disc only in 
those regions where the coolign rate A exceeds the adiabatic 
heating rate (see the discussion in DV09) . This condition de- 
termines the location of a radius Rsf where stars are allowed 
to form. 

We denote the time lag between cloud formation and 
collapse to form stars as tiag. This is basically the lifetime of 
molecular clouds; in our model it corresponds to the delay 
between the onset of instability (that we assume to be the 
driver leading to the formation of the clouds) and the start 
of the star formation event. Observational estimates for lag 
times range from 1 to 7 Myrs (Klessen et al. 2009). We 
assume tiag — 3.5 Myr. We have checked that our results 
are not strongly affected by the assumed value of fiag: final 



the Lyman- Werner (LW) band (11.2-13.6 eV). LW photons can 
photo-dissociate H2 molecules, thus suppressing cooling below 
8000 K, as dictated by Lyo cooling. Gas could in principle cool 
down to much lower temperature in presence of H2 - Less massive, 
colder discs form, leading to a much lower inflow. This prevents 
the formation of BH seeds unless higher inflow efiiciencies are 
considered. We will consider this more general case in our next 
paper, where LW flux ensuing from star formation history in a 
cosmological context will be calculated. 
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black hole seed masses change by a factor 2-3 for tiag ranging 
between 1 and 7 Myrs. 

We assume a star formation rate surface density Esfr 
that follows the empirical Schmidt-Kennicutt law (Kenni- 
cutt 1998): 



EsPR = 2.5 • 10" 



E(i?) 



IMqpc- 



Meyr ^kpc 



(7) 



and calculate the star formation rate J\f,,d by integrating 
over the region of the disc where stars form. We thus inte- 
grate SsFR between Rtr and Rsf to infer a star formation 
rate M, d in the outer disc: 



M*,d = 57r X 10"* [Eoiio]'''^^ (K 



,3/5 
'•SF 



3/5n 



(8) 



Star formation decreases the amount of gas that can be 
funnelled in the inner disc. For this reason and as long as 



> Merit we assume a net inflow rate Mini equal to 



Mint = Mgrav - M, 



(9) 



2.4 Supernova feedback 

Approximately 3 Myrs after an episode of star formation, 
the most massive stars explode as SNae. Metals processed in 
their centers are released into the surrounding gas, increas- 
ing the metal content. Metal yields and energy production 
depend on the initial mass function (IMF) of the stars. We 
here assume a Salpeter IMF with minimum (Mi) and maxi- 
mum (Mu) mass of 0.1 and 100 Mq respectively normalized 
to 1 Mq. 

We correlate the rate of metal production with the star 
formation rate in the disc assuming that the total mass in 
metals scales as: 



Mmct = YynctlyST^M^ 



(10) 



where Ymct is the total IMF-averaged metal yield and usn is 
the fractional number of SNae exploding after the formation 
of a mass M, in stars. We assume that SN progenitors are 
stars with masses between 10-50 Mq (Scannapieco et al. 
20030. With our prescriptions Kmet = 1.6 Mq and usn = 
0.00484. 

SN explosions in high redshift halos can be extremely 
destructive. SN driven bubbles can push away part of the 
gas in the halo, eventually fully depriving the host of its 
gas. The efficiency of gas depletion depends on the depth 
of the potential well of the host and on the energy released 
during the explosion. 

Using energy conservation we compute the amount of 
mass that is removed from the disc plus halo as a result of 
SN esplosions, assuming that each single explosion releases 
an energy _Bsn = 10^^ erg (Woosley & Weaver 1995). The 
resulting outflow rate M^h is 



^ We here follow Tsujimoto et al 1995 and assume that stars with 
masses greater than 50 Mq collapse into a BH without ejecting 



2 (IH- fib/firn) GMh/i?vir + ^V^, ' 



(11) 



where is the fraction of energy channelled into the out- 
flow, is the velocity of the outflow at the virial radius. 
We refer the reader to Appendix A for details of our recipes 
for gas depletion. Lastly, we define the retention fraction co- 



efficient fri 



(Mg 



Afsh)/Afgas. After the explosion we 



assume that a fraction /^ot of the metals yields produced 
remains in the host. 



2.5 Disc evolution 

The evolution of an unstable disc is followed tracing the time 
dependence of Mcoid, Mint, M*,d and Msh, integrating equa- 
tions 2, 8, 9 and 11 assuming that initially all the gaseous 
mass is hot (T — Tvir). 

Changes in these quantities over cosmic time correspond 
to changes in the scale parameters of the Mestel disc (Eq, 
_Ro) and _Rtr. Angular momentum conservation and total 
mass conservation in the system are imposed to infer the 
following equations linking the current values of Eo(t) and 
Ro{t) to their value at the moment in which the disc reaches 
the instability regime: 



Eo(t) = Eo(0) 



Roit) = i?o(0) 



1 



md{t) 



(12) 



(13) 



heavy metals. 



where md{t) — Mcooi{t)/Mh and ma{t) is the total frac- 
tional gaseous mass removed from the outer disc, either via 
star formation, inflow and outflows from SNae. The disc al- 
ternates episodes of instability to states of quiescence: as Q 
increases above Qc, the instability is quenched. Disc param- 
eters evolve regulated by the amount of new mass that is 
added by accretion from the halo. 

We halt the simulations as soon as the hot halo gas is 
consumed and the disc is in a stable state. We proceed on 
studying in the next section the evolution of the gas that 
has been funnelled at the center of the disc where we expect 
that a NC and/or a BH seed form. 



3 CENTRAL MASSIVE OBJECT FORMATION 

In the previous Section we focused on the evolution of 
Toomre unstable discs. Here we study the formation of a 
NC from a mass reservoir Mini, funnelled into Rtr after fiag. 

Star formation in the central part of the disc can pro- 
ceed with a higher efficiency with respect to the outer disc as 
higher pressure and density are expected there (Elmegreen 
& Efremov 1997, Li et al. 2007). We assume that a frac- 
tion esF of Minf is converted into stars and we here adopt 
esF = 1 (Li et al. 2007). The mass and radius of the cluster 
are Md = espAfinf and Rd = min(_RsF, -Rtr). 

NCs formed in this way are very compact and unstable 
objects. Numerical studies of stellar clusters have shown that 
gravitational encounters between stars can drive the clusters 
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towards core collapse, leading to a contraction of the core 
where physical collisions between stars are frequenlQ. Core 
collapse is driven by dynamical friction of the more massive 
stars. Star-star collisions between massive stars in the cluster 
core can happen in a runaway fashion, leading to the build 
up of a very massive star (VMS), whose mass can be as high 
as 1000 M0 (Portegies Zwart et al. 1999,2004; Freitag et al. 
2006a,b; Gurkan et al. 2004). 

After ~ 3 Myrs the first SNae explode. For the systems 
we are interested in, this happens before the inflow process 
ends. We assume that a NC is able to produce a BH seed only 
if it undergoes core collapse before the first SN explosion. 
We compute the mass of the nuclear star cluster, M^ifi, 3 
Myr after the onset of star formation, and from these we 
infer the core collapse timescale, tec- For those systems with 
tec ^ 3 Myr, we calculate the mass for the very massive star 
MvMS as in DV09: 



Ad- 



VMS 



: m. + 4 X 10"^Afci/c In Ac 



^^j^_3Myr^ 



tc 



(14) 



where — 0.2 and In Ac ~ log(0.1A/ci, 0/(^1)) (Portegies 
Zwart & McMillan (2002)). 

At low metallicities, stars more massive than 260 Mq 
are expected to leave a remnant BH, that retains most of its 
progenitor mass. Each time Mvms is greater than 260 Mq 
we assume a BH seed is formed with Mbh ~ A/vms. For 
a BH seed to form the stellar cluster needs a total mass ^ 
260 Mq in massive stars that can segregate. We here assume 
that stars more massive than 10 Mq efficiently segregate 
in the cluster core. Once all these more massive stars are 
consumed collisions are assumed to stop. For our chosen 
IMF a minimum VMS mass of 260 Mq, corresponds to a 
minimum cluster mass of 1.5 X 10* Mq. We here include 
this caveat and assume that BH seeds can form only in those 
nuclear clusters more massive than this threshold. 

The ability of runaway collisions to form an interme- 
diate mass BH also depends on the metallicity of the gas 
from which the stars originated. Recent simulations of star- 
star collisions have shown that at metallicities higher than 
10""^ Zq, physical collisions are too disruptive, and the build 
up of a star with mass greater than 260 Mq is inhibited 
(Glebbeek et al. 2009). For this reason we assume that a 
BH seed forms only at metallicities < 10~^ Zq. 

Summarising, BHs form provided that the following 
conditions are met: (i) Z > Zcrit so that a NC can form; 
(ii) A is sufficiently small for the disc to reach gravitational 
instability; (iii) clusters form with tec < 3 Myr and (iv) 
massive enough to host an adequate amount of stars with 
m* > 10 Mq. 



* Note tliat these studies concentrate on halo stellar clusters. NC 
dynamics can be alterated with respect to their halo contcrparts 
by (i) the presence of the external stellar environment and (ii) 
continuos gas accretion. These two effects work on opposite di- 
rections: cluster contraction can be reverted as stars of the more 
rarefied outer-system arc trapped into the central clusters. In- 
flow of gas can instead drive an acceleration of core collapse as it 
deepens the potential well. The final outcome resulting from the 
competition between them is not easy to infer and it is behind 
the aim of this paper. We will address this issue in a follow-up 
paper. 
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Figure 1. Upper panel: Mi^f (solid line) and M<, ^ (dash-dotted 
lines) versus time t for our reference models. The triangle high- 
lights the BH mass at the moment of its formation. Lower panel: 
Afinf (solid line) and A/, ^ (dash-dotted line) for our reference 
model as a function of time. The horizontal dotted line marks the 
critical threshold for fragmentation Merit ■ Vertical lines define the 
different evolutionary phases of the disc properties as described 
in the text. 



4 RESULTS 

In this Section we describe how disc properties are shaped by 
gravitational instabilities and star formation. Our reference 
halo has mass Afh = 6.5 x 10^ Mq, and forms at Zvir = 10 : 
we further adopt A = 0.05, Z/Zq = 10"", 77 = 0.3, Qc = 2, 
esF = 1 and tiag = 3.5 Myr. 
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4.1 Disc and cluster properties 



In this subsection we illustrate the evolution of the reference 
disc model. After a time rinst = 1.3 x 10** yr, the disc has 
grown enought to hit the instability threshold as it has as- 
sembled enough mass so that the Toomre parameter drops 
below Qc- The upper (lower) panel in Fig. 1 shows Minf 
(Minf) and M«,d (Aft,d ) as a function of time t, from the 
time the disc becomes unstable until the time of formation 
of the central BH. The evolution proceeds along different 
phases: 

1) Tinst < t < Tinst + tlag: at f = Tinst the disc starts col- 
lecting gas into the center and this leads to an enhancement 
of Minf. The disc parameters are such that Mgrav is greater 
than Merit (horizontal dotted line in the righ panel of the 
Figure) . Clump formation in the disc sets in but these newly 
formed clouds have not had enough time to collapse until a 
time tiag has elapsed. 

2) Tinst + tiag < t < Tinst + tlag + 3 Myr: Star formation be- 
gins. The critical density for star formation remains constant 
as metals have not been released yet and the star formation 
radius still remains close to its initial value. Both Afinf and 
M«^d increase until the first SNae explode. We label the nu- 
clear mass at this time as M^ifi. 

3) t > Tinst + tiag + 3 Myr: the first SNae explode affecting 
the subsequent evolution of the system. Part of the metals 
produced in stars are released in the surrounding gas in- 
creasing the metallicity above Za-it- At the same time SN 
explosions evacuate gas from the system lowering the sur- 
face density. As a result the Toomre parameters increases 
above Qc. 

The disc then lingers on a marginally stable configura- 
tion thanks to the high Msh that is able to counter-balance 
Mcooi- Fragmentation in the disc is temporary halted as 
Mgrav drops below Merit. At later times, the surface den- 
sity starts to increase again, as hot gas in the halo continues 
to cool. The disc then enter a new cycle of instability, and 
evolution terminates when the hot gas reservoir has been 
consumed. 



4.2 Dependence on input parameters 

We here explore disc evolution as a function of the halo 
mass, Mh, for three different values of inflow efficiency, 77, 
and as a function of the initial metallicity Z. 

Fig. 2 shows the NC mass Mci,o (solid lines) and BH 
mass Mbh (log dashed lines) at the time of BH formation, 
and the stellar disc mass M«,d (dash-dotted lines), and net 
inflow mass Mint (dotted lines) at the end of the simulation, 
as a function of i\fh, for different values of rj. 

The NC mass, Mci,o, depends weakly on the host mass, 
Mh, and somewhat more strongly on the inflow efficiency, 
77. By contrast, the final stellar disc mass depends signifi- 
cantly on Mh. This is because (i) more massive halos have a 
larger gas mass and retention efflciency (/ret); (ii) gas cools 
more rapidly in larger halos. This allows for a larger depar- 
ture of the Toomre parameter Q from Qc, and larger inflow 
(Mgrav, see equation [5]). Mgrav is larger than the critical in- 
flow triggering fragmentation. Merit, and more stars form in 
the outer disc. For 77 = 0.1 a threshold in Mh exists below 
which Mgrav < Merit always, and A/*,d drops to zero. 
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Figure 2. Mass M^ifi (solid lines) of the NC and Mbh (dashed 
lines) at the time the BH forms, and Minf (dotted lines), M, d 
(dash-dotted lines) at the end of the simulations, as a function of 
A/h for different inflow efficiencies. High, central and low panels 
refer to r; = 0.1, 0.3 and 1, respectively. Reference properties for 
the halo are 2:vir = 10 and A = 0.05; for the baryonic component 
Z/Zq = 10"*, Qe = 2 and esp = 1. 
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Figure 3. M^i q (solid linos), JV/bh (dashed lines), Mjnf (dotted 
lines) and ^ (dash-dotted lines) as a function of the initial 
metallicity of the gas for our reference halo. The halo mass is 
kept constant at its reference value A/h = 6.7 X 10^ Mq. Vertical 
lines mark the minimum metallicity for low-mass star formation, 
Zcrit — IQ-^'*'^ Zq , and the maximum metallicity that allows 
BH formation, Z^rit — lO"'^ Zq. 

Dotted lines in Fig. 2 show the net inflow mass, Minf, 
at the end of the simulations. Mini increases with increasing 
Mh as a result of both initially larger gas masses and deeper 
potential wells. Minf depends weakly on rj. For higher rj more 
gas is consumed by extended star formation before it can 
reach the nucleus and final Minf are lower. 

Disc evolution depends also on the initial metallicity Z. 
The metallicity determines when stars can form and it is 
thus a key parameter: there is a minimum metallicity where 
'normal' stars star forming, Zait ~ 10-4-87 and a maxi- 
mum metallicity that allows BH formation, Zcth — Zq. 

In Fig. 3, the NC mass Mci_o (solid line), BH mass Mbh 
(dashed line), net inflow mass Minf (dotted line), and stellar 
disc mass M*,d (dash-dotted line) are plotted as a function 
of the initial gas metallicity. Star formation sets in only for 
Z > Zcrit (left vertical dashed line in Fig. 3). For the major- 
ity of the systems we are considering here, Rsf > Rti- The 
flrst three quantities are almost independent of metallicity, 
as long as Z is in the range between the minimum and max- 
imum value. The disc stellar mass is instead very sensitive 
to the exact value of the metallicity. A'h^d increases with Z 
(see equation [8|. Mass is consumed instead of flowing into 
the nucleus, and Mci,o weakly decrease^. When the metal- 
licity is large enough that Rsf equals Ro than the whole disc 
forms stars, and the dependence on Z disappears. 

In our reference model we set A = 0.05. The distribu- 
tion of spin parameters can be described as a log-normal 
distribution with peak at A ~ 0.03 — 0.05 and dispersion 
~ 0.5. For our reference halo the maximum A allowing for 



instability is 0.07. Approximately 96% of the halos in the in- 
teresting mass range have A < Amax and are therefore prone 
to disc instabilities. 

If Qc is lowered to unity, instabilities are easily 
quenched, and the process becomes less efficient (see upper 
right panel in Fig. 4). 

4.3 NC and BH masses versus redshift 

In this Section we apply our prescriptions for the cooling of 
gas, star formation and inflow to determine the properties 
of the NC and BH populations, in low metallicity environ- 
ments. 

To this purpose, we construct a sample of halos with 
the same properties of those analysed in DV09. We focus 
on halos with Zvir ranging between 5 and 30 corresponding 
to Tvir between 10* K and 1.8 x lO" K (Barkana & Loeb 
2001), and we determine their frequency using a modifled 
version of the Press & Schechter formalism (Sheth & Tor- 
men 1999) in a WMAP5 cosmology (Dunkley et al. 2009). 
To each halo, we assign a value of the spin parameter. A, 
extracted from the probability distribution found in the Mil- 
lennium simulations. We here adopt the same prescriptions 
for the metallicity evolution as in model A in DV09, while 
we adopt the model described in this paper for the evolution 
of each disc-halo. The model is thus an improvement of that 
discussed in DV09. More realistic models that follow metal 
enrichment self-consistently will be described in Paper II. 
We consider different inflow efficiencies »7=1 (dotted lines), 
0.3 (solid lines) and 0.1 (short-dashed lines) for Qc ~ 2 and 
1] — 0.3 for Qc~l (long-dashed lines). 

In the upper left panel of Fig. 4 we show the evolution 
of the NC mass, {Afci,o) (averaged over the sample) as a 
function of redshift. NCs start to form at z = 14, as halos 
are polluted above the critical metallicity for fragmentation 
(see DV09B BH seeds start to form after the formation of 
the flrst NCs. In the upper right panel of Fig. 4 we com- 
pute the fraction /bh of halos with Tvir > W K hosting a 
BH seed, as a function of redshift. /bh is regulated by the 
fraction of halos with Zcrit < Z < IQ-^ Zq: it first mcreases 
as more and more halos are polluted, and then decreases 
when most halos reach Z = W~'^Zq. /bh in this paper are 
larger than those found in DV09 (see their Figure 8, upper 
panel) by a factor 10. With the prescriptions we adopt here, 
Mcooi evolves with time, regulated by cooling of the halo 
gas whereas in DV09 a fixed Mcooi/Mh is used. This implies 
that for a given Qc, the fraction of unstable discs that form 
is increased. 

Mean black hole masses (Mbu) versus redshift z axe 
plotted in Fig. 4 (lower left panel). (A/bh) spans values be- 
tween a few hundred Mq and 2000 Mq. For most of the 
NCs considered here, the growth of the VMSs is limited by 
the amount of stars with m* > 10 Mq. Lower (Mdfi) thus 
lead to lower mass seed BHs. An exception occurs for the 



^ Note that for systems with ijgp < iitr, the initial cluster mass, 
A'/(,i 0, increases with Z as (i) extended star formation is not al- 
lowed in the outer disc, so that the amount of mass funnelled into 
the nucleus remains constant; (ii) higher mctallicities correspond 
to higher -Rsf so that to a larger amount of nuclear gas can be 
converted into stars. This case was considered in DV09, we refer 
to that work for a discussion of this regime. 



^ In our reference model (r? = 0.3, Qc = 2) the mean cluster 
mass {Mcifi) decreases with redshift as the metallicity increases. 
This behaviour appears only as the majority of the systems we 
are considering here have i?SF > ^tr (that is, stars form in most 
of the disc). If this were not the case, the behaviour of (Af^i o) 
versus redshift would be the opposite (see DV09). 
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Figure 4. Panels from up to low: (A/d o) (for those clusters able to form a seed BH), /bh> (A/bh) and Psccd a-s a function of rcdshift 
for different values of rj and Qc- Line styles: solid lines refer to a model with rj = 0.3, Qc = 2, dotted lines refer to a model with rj = 1, 
Qc = 2, short-dashed lines refer to a model with rj = 0.1, Qc = 2, long-dashed lines refer to a model with r] = 0.3, Qc = 1. 



highest Mcifi (i.e. 77 = 1): for these systems tec is long and 
not all massive stars contribute in building up the VMS. 
(Mbh) stays at a near constant value ~ 900 Mq, at any 
redshift. 

Comoving BH seed densities are plotted in the lower 
right panel, pscod depends on both BH formation efficiencies 
and masses attained. They are found to depend weakly on 
the parameters: final Pscod only differs of a factor 2, ranging 
between 2000 and 4000 M© Mpc"^. 



5 SUMMARY 

In this paper, the first of a series, we presented a model for 
the formation of high redshift, metal poor NCs and their 
embedded BH seeds. The channel that we propose can be 
summarized as follows: 

• After the virialisation of the halo, hot gas cools down 
and the presence of angular momentum leads to the forma- 
tion of a cold gaseous disc. The disc accretes gas from the 
halo and grows in mass, increasing its surface density to the 
extent that gravitational instabilities develop. This leads to 
the growth of non-axisymmetric perturbations that redis- 
tribute angular momentum, driving mass into the central 
part of the halo. 



• For those systems where the metallicity Z exceedes Zcrit 
and Mgrav > Merit, stat formation in the extended disc can 
set in, consuming part of the mass that would otherwise been 
accreted in the center. The inner density profile steepens 
inside a transition radius where a compact nuclear stellar 
cluster forms. 

• In those NCs where tec < 3 Myr, core collapse antic- 
ipates stellar evolution and runaway collisions between the 
most massive stars build up a very massive star that finally 
implodes leaving a remnant BH seed. 

• The evolution continues and it is driven by the competi- 
tion between star formation, triggered by prolonged cooling 
of hot halo gas and SN explosions that evacuate part of it. 

• Evolution halts as soon as the hot halo gas is con- 
sumed/blown away and the disc settles into a state of sta- 
tionary equilibrium. 

We apply our recipes to a set of halos with Tvir > 10* 
and properties typical of z ~10-20. We determine how 
changing halo properties and model parameters influence 
the final NC and BH seed population. The NCs have masses, 
at the moment of BH formation, of 10* — 10^ Mq and host 
BHs with masses between 300 and a few 10^ M©. 

We compared our new prescriptions with the results of 
DV09. The mean NC masses are found to decrease with 
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decreasing redshift: higher metaUicities usually correspond 
to larger star formation rates in the extended disc. These 
reduce the amount of gas channelled into the center, leading 
to less massive clusters. Except for the more massive NCs, 
Mbh is dictated by the availability of massive stars. The 
decrease in Mci,o consequently leads to less massive black 
holes. The total mass densities in BH seeds, Pseed, are higher 
with respect to DV09 as a result of the larger number of 
gravitationally unstable discs. In our model we do not take 
into account any BH growth by accretion, nor their assembly 
during galaxy clustering. Comoving densities found here are 
thus lower limits. 

In the second paper of this series we will apply this semi- 
analytical scheme to merger tree histories to trace NC and 
BH seed formation into a cosmological context. In the third, 
we will follow the evolution of the emerging BH population 
through accretion and clustering and compare our results 
with current observation of BH occupation frequency in low 
redshift galaxies, quasar luminosity functions, etc. 



APPENDIX A: GAS DEPLETION FROM SN 
FEEDBACK 

In this section we describe the prescription used to infer the 
amount of mass lost from the host following a SN episode. 
Under the assumption that the gas follows an isothermal 
profile 



Pgas — 



(Al) 



we compute the initial binding energy of the system Sb.in 
inclusive of the dark matter and gaseous componenlQ. 



p _ GMl ( fib 

iib,in ^ 1 1 + TT" 



(A2) 



After removing a gaseous mass Msh the new binding 
E^t is 



a mass M, of stars, divided by M«. SNae start to explode 
after 3 Myr from the onset of the star formation episode. 
We assume that at that time all stars with masses in the SN 
regime explode together giving rise to a single bubble and 
that a fraction of the total energy released by SN explo- 
sions, is channeled into the outflow. Following Scannapieco 
et al. (2003), we relate with the halo mass as follows. We 
define (5b (Mh) as 



1.0, 



5b(M) = <^ 1.0 - 0.165 In(iV-i), 

\ [1.0 - 0.165 In(lOO)] imN^^ 



Nt s; 1 

1 < iVt ^ 100 

100 s; iVt 

(A5) 

where iVt = 1.7 x 10"^(f7b/t^m)Mh/ M©. /„ is then given by 
f-u, = 5B(Mh)/&(Mh = 2x 10* Mq). This prescription fol- 
lows from the scaling relation found in Ferrara et al. (2000) , 
normalised with the result found by Mori et al. (2002). 

Let us now compute ifsN = l/2MshVsij. Here the shell 
velocity is calculated at the moment it reaches -Rvir. The 
topology of early metal enrichment can be extremely com- 
plex. Enriched gas first propagates into the cavities created 
by the shock (Greif et al. 2007), and after the bubble shell 
leaves the host, it spreads mainly into voids, leading to a 
very anisotropic spatial distribution of metals. We here ne- 
glect these degrees of complexity and simply assume that 
SN-driven bubbles evolve in spherical symmetry. The shell 
radius is assumed to follow the Sedov- Taylor solution, as 
long as the pressure of the bubble drops below the value of 
the surrounding medium. The shell radius R^h evolves with 
time t as 



Rsh = 1.15 



Pbk 



1/5 



(A6) 



where t is the time after the explosion and pbk is the am- 
bient gas density, scaling as (r2b/f2m)Afh/-Rvir- The velocity 
of the shell v^h in the Sedov- Taylor approximation is given 
by 2/5-Rsh/t. Under these conditions, Vsh{Rvir) can be easily 
inferred. At the time R^h ~ Rvii the kinetic energy is 



b,f 



Rvii 



1 + ^-^) • (A3) 



Energy conservation imposes the following equality: 
Eh,in + £sN = -Eb,f -I- KsN, where Ssn is the energy injected 
by the SN explosions and Ksn is the kinetic energy of the 
outfiow that is leaving the halo . 

The energy £sn by SNae, over a time t, is given by 



£sN{t) = I'sN-EsN / M«ds 
Jo 



(A4) 



where i?sN is the energy of a single SN explosion, set equal to 
10^^ erg, M* = espMinf -I- Mt,d and usn is the ratio between 
the total number of SNae exploding after the formation of 



^ Note that the actual gas density profile at the moment of SN 
explosions would be quite different from the initial one. Difl'cr- 
ences would arise both from the fact that the gas needs to cool 
down in order to form stars, and from the radiative feedback of 
stars during their lifetimes. 



AsN = 0.67-— ——/„isN. 

"b 



(A7) 



To compute A4sh for a given halo we recure to energy 
conservation leading to 



1-0.67 



Msh \ 



fib MkJ 



fwSsN ~ 



Rvir 



2 1 



Hb \ _ M,h " 

fimy Mh 

(A8) 
less than 



We note that Mgh/Mh is always much 
2 {1 + Qb/flm)- Negletting the term Mgh/Mh in parentesis 
on the right side of the equation the outflow mass of the 
first single bubble is given by 



M,h 



fm£sN 



(A9) 



2 (l-f Ob/finO G Ml, /Rvir + 

We assume that an outflow develops only when v^h > Vcsc- 
Otherwise Mah ~ 0. 

Massive halos are not compleately evacuated after this 
first SN episode. Further star formation is often allowed and 
to account for subsequent SN esplosion episodes we compute 
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the total number of exploding SNae A^'sn between time t\ 
and t2 



«^SN UspMinf + Af,.d ds 



(AlO) 



We define Msh,s as the contribution of a single esplosion 
(i.e. equation IA9I with £sn = ^^sn). An Outflow develops 
only if it injects sufficient energy so that > Vcsc- This re- 
quires a minimum energy in a single bubble, i.e. a minimum 
number of SNae. For a given M* , we calculate the expected 
number of SNae on a timescale of the order of their progen- 
itor lifetimes (i.e. 3 Myr). If the resulting v^h is greater than 
the escape velocity from the halo, outflows are allowed and 



diVsN 
' At 



(All) 



Otherwise M^h — 0. 
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